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Abstract 

Gravitational waves in the hnear approximation propagate in the Schwarz- 
schild spacetime similarly as electromagnetic waves. A fraction of the radia- 
tion scatters off the curvature of the geometry. The energy of the backscat- 
tered part of an initially outgoing pulse of the quadrupole gravitational radi- 
ation is estimated by compact formulas depending on the initial energy, the 
Schwarzschild radius, and the location and width of the pulse. The backscat- 
ter becomes negligible in the short wavelength regime. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Backscattering has been investigated for a long time for various wave equations (see, 
for instance, [jl|). In general relativity, this topic has been studied since early 1960's ( 0, 
0). This paper continues the programme that started with the study of the backscatter of 
scalar |^ and electromagnetic fields ( ||^ and [^]). Here we investigate the propagation of 
even-parity gravitational waves in a (fixed) background Schwarzschild spacetime, assuming a 
nonstationary source. The discussion, however, is done without any reference to the source. 
We only deal with field quantities. It is assumed that initial data are those of an isolated 
pulse (burst) of a gravitational wave. The main question that is answered is what fraction 
of the initially outgoing radiation may undergo backscattering before reaching null infinity? 
The strength of the backscattering is assessed by bounding the fraction of the initial burst 
energy that will not reach a distant observer in the main pulse. 

The even-parity waves are the only waves which are radiated during the axisymmetric 
collision of nonspinning black holes and since in this case the Schwarzschild spacetime 
is a valid starting point for an approximation scheme, it gives us an opportunity to bound 
the strength of the phenomenon in a fairly realistic astrophysical context. 

The following five sections of the paper give a theoretical description of the backscattering 
effect. The Sec. II brings notation and the Zerilli equation. In Sec. Ill, the initial data are 
bounded by the initial energy and the solution is sought in the form of superposition of an 
outgoing radiation (defined by initial data) and a backscattered term. The evolution of the 
backscattered term can be bounded by solutions of two differential inequalities. The bounds 
that are derived in Sec. IV and in the Appendix deal with a general situation; no assumption 
is made about the initial radiation. In Sec. V, we discuss initial data that are of compact 
support and in addition, the relative width of the support is small. Such data correspond 
to radiation that is dominated by short wavelengths. In this case stronger estimates are 
derived. They imply that in the limit of short wavelengths (relative width of the support 
tending to zero) the backscattering effect becomes negligible. In Sec. VI, the "small relative 
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width condition" of Sec. V is supplemented by the assumption that the initial burst is far 
away from the horizon. 



II. FORMALISM 

The space-time geometry is defined by a Schwarzschildean line element, 

ds^ = _(i _ ^)rft2 + -A^dR^ + R'dn' , (2.1) 

where t is a time coordinate, is a radial coordinate that coincides with the areal radius 
and dQ"^ = dO"^ + sin^ Odcj)"^ is the line element on the unit sphere, < < 27r and < < vr. 
Throughout this paper G, the Newtonian gravitational constant, and c, the velocity of light 
are put equal to 1. 

As explained in the Introdction, we restrict ourselves, to the even-parity axial pertur- 
bations. Their propagation is ruled, in the linear approximation, by the Zerilli equation 
0. Formulated in terms of the gauge- invariant amplitude defined by Moncrief |p, this 
equation reads, in the case of / = 2 multipole M, 



-d^ + d^,,)^ = V^, (2.2) 



where the potential V is given by 



V(m - fin 2^^^ ^ + n 2m 63m^(l + f) 



and where 



r* = i? + 2mln(— - 1) (2.4) 
2m 

is the tortoise radial coordinate. 

Consider a set of functions "if i{r* — t), i = 0,1, 2, that satisfy the following linear relations 

dt^i = 3^0, 

dt^2 = *i - mSt^i. (2.5) 



The combination 

^.Mr'-t).^^.^^ (2.6) 

solves Eq. ( |2.2| ) in Minkowski space-time (m = 0); it represents purely outgoing radiation. 

Let the initial data of a solution of ( p.2| ) coincide with ^ at t = 0. Then, initially, \E' 
represents a purely outgoing wave. It should be noted that the assumption that initial data 
are (initially) purely outgoing is made in this paper only for the sake of clear presentation. In 
the linear approximation the propagation of the initially outgoing radiation is independent 
of whether or not ingoing radiation is present. 

We decompose the sought solution \E'(r*,t) into the known part and an unknown 
function 6 

^ = ^ + 6. (2.7) 
Due to the choice of the initial data made above one has 6 = dt6 = 0, at t = 0. 



III. ENERGY ESTIMATES 

Let us assume that the quadrupole initial data are defined by a smooth triad of the 
functions \E'fc {k = 0, 1,2) with the initial support [a,b] (h < oo). That guarantees that the 
initial energy density multiplied by i?^, p = {{dt'^Y + {dr*"^)^ + V'^^)/riji, is smooth and 
vanishes on the boundary a. Here rjR = 1 — 2771/ R holds. 

The energy content inside a part of a Cauchy hypersurface St that is exterior to a ball 
of a radius R can be defined as E{R,t) = J"^ drp{r,t). Let us point out that in order to 
ensure a proper normalization of the energy flux at infinity, there should be a normalization 
constant in the definition of the energy E{R, t). We decided to omit it, since later on we shall 
be interested only in the relative efficiency of the backscatter and thus the normalization 
factor cancels out. The total initial energy corresponding to the hitherto defined initial data 
is equal to E{a, 0). 
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Lemma 1. Defining 



Ci = |j(2 + y273)£;(a,0), 



2 



= V2^(a,0), 



2E{a,0) 



and introducing the two nonnegative functions 

, f-2m + R\ J -1 1 \ 



(-2m + i?)4 (-2m + a)V V3(-2m + i?)3 3{-2m + a)\ 

20m2( — — ^— - + -— i 1 + 10m ~^ ^ + ^ 



(-2m + i?)2 (-2m + a)2y V-2m + i? -2m + 



and 



4/ -1 1 



^^^""'^ ((-2m + i?)2 + (-2m + 0)2) + 

(24m^ + +8mln 

^ ^V-2m + i? -2m + aJ \-2m + a^ 

the following inequahties hold at i = and for it! > a: 

^ < c.,f ^m. 



\J a a 



i? - V 'R 



Proof. 

One can explicitly verify that 

The equations (^.5|), using the relations ([^■6D, result in 

^ / -3/2 ^1 \ 3 / dt^ ~ 2^y^\ 

^ / -1/2^^2 X 1 /^gt^ , , 3m ^ 

The integration from a to R and the use of the Schwarz inequality yields 

Integrating J^dr:^, one immediately arrives at the first of the postulated inequalities. 

In order to show the third inequality, notice that |^(-R)|-R^ = | drdr{'^{r)r^)\. The 
latter expression is bounded from above, using the Schwarz inequality, by 



ma,0) RVE.oJ,_f-f-A'^\ (3.8) 



where the inequality in Eq. (|3.8|) follows from the monotonicity of the energy as function 



of R. The first factor on the left hand side of this inequality is not greater than y2E(a, 0) 
since Qrir/r"^ < V{r)/r]r. The replacement of r/^^ by 7]^^ and the integration of the other 
factor leads to the desired result. 

The second of the equations of (|3l6|) can be integrated. The Schwarz inequality and 
direct integration as well as the \I' and \E'i estimates should be used in order to get the 
second inequality of Lemma 1. The \E'o-estimate, in turn, follows from the identity \E'o = 
\E' — \I'i/i? — '^2/R^ and the preceding estimates. 
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IV. THE ESTIMATE OF THE DIFFUSED ENERGY 



Let us define the strength of the backscattered radiation that is directed inward by 



h^R, t) = —{dt + dr*)5{R, t). 

m 



(4.1) 



Let the outgoing null geodesic f(R,t) originate at {R,t). If a point hes on the initial hyper- 
surface, then we will write T{Rfi) = Lr. By r(7^Q (^ t) will be understood a segment of 
f(Bo,to) ending at {R,t). 

A straightforward calculation shows that the rate of the energy change along is given 

by 



{dt + dr*)E{R,t) = 
- rilh'_{R,t) + VS\a) 



(4.2) 



It is necessary to point out that in the case of the initial point Rq > a the result would be 
more comphcated; the differentiation of the energy along would depend also on 
and "^2- If, however, the outgoing null geodesies is r^, then it starts from a where \E'o, ^E'l 
and ^2 do vanish. Since these functions depend on the difference r* — t, their values along 
outgoing geodesies are constant, and that allows one to conclude that they vanish at 

The energy loss, that is the amount of energy that diffused inward is equal to a line 
integral along r^. 



SE. 



f 

J a 



a 
oo 



dr 



E{a, 0)-E^^ 
V5''' 



r]rh_ + 



(4.3) 



Our goal is to find an estimate of 5Ea of a single pulse of radiation basing only on the 
information about the position and the energy of the initial pulse. Obviously, < 5Ea < 
E{a, 0) holds. We are interested in deriving in this section a frequency-independent bound, 
but later we obtain estimates that are frequency-sensitive. 

S is initially zero and its evolution is governed by the following equation 



-d^ + dl,)5 = V5 + iV- 6^) (^0 + § + §) + 



\I/2 

-3^1 + 2-^ 



(4.4) 



One can define an "energy" H[R,t) of tlie field 6 which is contained in the exterior of a 
sphere of radius R as follows 



H{R,t) 



h r]r{drO) + — 

Vr 



The rate of change of H along T(^R^t) is given by 



(4.5) 



{dt + d:)H{R,t) 



-Vr 



Am 



V 



Vr 



R 



d,6 



Vr 

2^2 63m(l + ^)/ ^1 



4(l + fi)2 



4m / dr^— 

R r'^ 



2^2 63m(l + 
-3^1 + ^+ ^ 



2r 

m ' 



4(1 + ^)2 



^ ^1 ^2' 

^0 + — + 



< 



(4.6) 



Herein the inequality follows from the omission of the nonpositive boundary term. This 
allows one to estimate the maximal value Hm of the 5-energy if, namely 



— mJE(a,0) 
Hm < 10.43 V ^ ' ^ ^ Q(^^2y 



(4.7) 



The calculation is essentially simple, but the algebra is quite lengthy and some numerical 
integrations are required. Details are relegated to the Appendix. We would like to point 
out that the 0{m?) terms become dominant only when the location of the initial radiation 
pulse is smaller than 6.6m. At a = 15m the neglected terms contribute much less than the 
leading term proportional to m. 

Now, the integration of the first part of Eq. ( [4. 61 ) along T(^a,Q) yields 



dR 
dRAm 



Vr[— + dR5 
^Vr 



a 

00 roo 



dr — - 



R 



-3^1 + 



Vr 
2^2 



63m(l + ^) 
4(1 + ^)2 



+ — + ^ 



r 



(4.8) 



Initially, H vanishes (both 6 and dt5 vanish) and H is manifestly nonnegative. The first 
integral on the right hand side of ([4 .81) is recognized to be just 5Ea. The second integral in 



turn can be shown to be bounded - using the Schwarz inequality and then the results of the 
Appendix - by 2^/^^(10.43 ""^f^"'"^ + 0{m^)). Thus, fj) implies 



8Ea < 2(10.43 V ^ ' ^ ^ 0{m^))^HM < 



2 



54.5 ( — 1 +0{m-^) 



E(a,0); (4.9) 



the right hand side of the first inequality achieves a maximal value when H is maximal and 
that implies the second inequality. 

Thus in summary, for the fraction of the energy that could diffuse through the null cone 
Ca, it holds: 

Theorem. 6Ea/ E{a,0) satisfies the inequality 

J^<54.5x(?^)%0(™V). (4.10) 

We would like to point out that the above derivation is more efficient and simpler than the 
one used in or when 6Ea was estimated directly on the basis of the estimates of 6 and 
h^. This alternative approach would require a laborious integration of the field equation 
and the final estimate would be much worse. 



V. THE WAVELENGTH OF THE INITIAL RADIATION AND THE 

BACKSCATTER. 

In this section we shall consider the backscatter of the radiation that is initially of 
compact support and, in addition, the condition {a — b)/a << 1 is satisfied. The leading 
contribution - only terms that are quadratic in - will be found. 

Under the above conditions one infers from Eq. ( |3.4| ) that on the initial hypersurface 

\^,iR)\0' < C,h'l'4^) < Ci^l^b'/' (5.1) 

and 
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\^2{R)\<C2b^Vb^ (5.2) 
are valid. With the same accuracy the inequahty ( p.3|) of the Appendix reads 



QmCnl^^y"'^ { dr^] + 4mC2b'Vb^( dr^] (5.3) 



^ ^ y^R J \Jr r^°y 

Herein, the integration extends from R, where R G r^, to R{b), which is defined by {R{b),t) G 
f b- One has R{b) — R = b — a up to the term m^. The integral /^^''^ dr^ is bounded from 
above by {b — a)/R^ and the integral fl^^'^^ (i^'p^ is bounded from above by {b — a)/R^^, again 
to lowest order in powers of m. 
Thus one arrives at 



(5.4) 



{d, + dr. ) ./mJ) < 4m{b - a)b^'^ + ^) . 

The integration of this inequality along the null geodesic Va yields 

^Hm < m{b - a)^(2Ci + C2) + 0{m'). (5.5) 
Taking into account the condition that b — a << a, one arrives at 

Since the amount of backscattered energy 6Ea is bounded from above by 2Hm, as shown in 
the preceding section, one finally arrives at the following estimate 

f b — a\'^ ( b\^ ( ?> 



5E^ 


< 


E{a,0) 




w? ( 


-) 




a I 



X - 

a J \a 

.3 I 3 



(2 1 ^ = 



2f ^^3 + 75 



+ 0(m7a'). (5.7) 



If {b — a)/a < 0.1, then the above formula predicts 



E[a,0) 
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It is clear that if the relative width of the initial pulse tends to zero then the effect becomes 
negligible. This can be translated into the dependence on the wavelength of the radiation 
0]: The compression of the support of a function leads to the decrease of its wavelength 
scale in its Fourier transform. 

A careful analysis of the higher order terms would show that they give a contribution to 
( ^.7|) that also scales like ^^^^ • '-^^^^ when a ~ 2m, Eq. ( |5.7|) would be of the form 

"%<CiM'-^t (5.9) 



E{a,0) - ^ '\2m 

where C{x) is a large number and x = 2m/ a. One can show that limj._,iC(x) = oo, but 
on the other hand C{x) is fixed, when 2m/ a is fixed. Thus ( ^.9|) implies that when b ^ a, 
then the backscatter becomes negligible. Radiation that is dominated by infinitely short 
wavelengths does not backscatter. 



VI. MORE ESTIMATES ON HIGH FREQUENCY RADIATION 



We assume initial data of compact support [a,b]. The initial energy E{a,0) (see the 
beginning of Sec. Ill) reads, expressed in terms of functions \l/o, "^i and \I^2, as follows 



no no 

E{a, 0) = / drp = dr 

J a J a 



-2\I^2(?^) - r(\I^i(r) + r(r\l/o(r) + \l/'i(r)) + \l/2(r)) 



(6.i; 



The radiation energy in the wave zone is known to be E{a, 0) = C /^^ (ir(\l/Q)^. This can be 
compatible with ( |6.1| ) (modulo a normalization constant, which is not relevant here), if the 
terms with (\I^o)^ give a leading contribution. 

One notices, that if "^^{R) (/i = 0,1,2) are of compact support, then |\l/^(i?)| = 



drdr"^ ^{r)\ < J{R — a) drdr^'iir). Combining this with (|2.5|) one arrives at 
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3|*o(r)| = |^'i(r)| < sVb^J J' dr{%) 



|^i(r)| = mr)\ < 2(6 - a)^/^ dri%y, 

\^2{r)\ < 0.8(6 - af'^slj\r{%Y. (6.2) 
Taking into account (p. 21), one concludes that if 

Ci(6-a)/a«l, (6.3) 
(Ci is a constant of the order of 100) then the energy is well approximated by -E(a, 0) = 

In such circumstances it is clear that our analysis can be greatly simplified. First of 
all the contribution coming from ^2 to the backscatter is much smaller than that due to 
\E'i; notice an additional power of (6 — a) /a in the relevant estimate of ( |6.2| ). Secondly, the 
estimate ( ]5.1| ) of \l'i is now replaced by a stronger result 

l^o(r)|<-^(6-a)3/V^(«'0)- (6-4) 

The repetition of the calculation of Sec. V gives finally (taking into account the above 
conditions) 

6E^ /2»w^y^ 



E[a,Q) \ a ) \ a 



VII. CONCLUSIONS 

In our paper we derived upper bounds for the backscattering of gravitational quadrupole 
waves propagating outward from a central compact object. The calculations were restricted 
to situations where the initial configuration was either an arbitrarily shaped wave with 
support outside some radius a or the wave was a sharp pulse, i. e. its extension was small 
compared to its initial location a. The obtained upper bounds show that, for a given central 
object, the backscattering is the weaker the more outside from the central object the waves 
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start propagating, and that is also the weaker the more compact the pulses are, i. e. the 
higher the involved frequencies are. The both results do confirm previous completely different 



calculations by Price, PuUin and Kundu |T^. Backscattering should thus be strongest for 
pulses which start propagating outward close to the horizon of a black hole. This claim, 
however, needs further investigation for the following reasons. First, we gave bounds from 
above and not from below for the amount of backscattered energy. Second, the linear 
approximation may not be accurate enough very close to the horizon. 

Results of Flanagan and Hughes fll]] and Buonanno and Damour have shown that 
the merger part of the gravitational wave signal could be a significant part of the total energy 
emitted. The wave pulse during the merger phase can be inside 3m. For a very compact 
pulse located in this region the inequality ( ^.9|) of Sec. V can still yield a nontrivial bound, 
but in the general case our estimates fail. The main reason why we are loosing much in the 
accuracy is that we are forced to use several times - for the sake of generality - the Schwarz 
inequality. The present bounds can be significantly improved if initial data are explicitly 
known, since in this case they can be numerically bounded by an exact expression involving 
the initial energy and the Schwarz inequality would be used only once. On the other hand, 
it has been discovered that the backscattering can be quite strong when a signal propagates 



from within the photon sphere [|I3 



In a forthcoming paper we shall discuss, and compare with the results of our present 
paper, several aspects of the backscattering of gravitational waves where the sources of the 
gravitational waves are taken into account. 
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VIII. APPENDIX 

In order to show the estimate ( [4.7| ) one begins with the second inequality of ( [4 .61 ). Notice 
that H{t = 0) = 0, since 6{R,t = 0) = dtS{R,t)\t=o- Thence the integration of ( |4.6|) along 
f a,(ij,t) yields 

1"^ dr , 

Vr 

where RHS{r) stands for the right hand side of Eq. ([4 .61) . Our task consists in estimating 
the line integral of \RHS{r)\. 

In order to do this one uses the estimates of (|3.4]). The calculation is quite long and we 
will describe only the main points. In the first step one uses the Schwarz inequality in the 
right hand side of (|4.6|), in order to obtain an expression of the type 



Hm< / -\RHS{r)\, (8.1) 



oo 



where /(r) denotes (— 3\E'i)^, (2\I'2/t)^, or the squares of the terms that are proportional to 



63m. The first integral can be bounded by y H{R), therefore (|4.6| ) and ( ^.2|) yield 

' .63m(l + ^)s2^2y/2 / 63m(l + ^^)s2^2\ 1/2 



The integrands of ( |8.3| ) are taken at a time t and {R, t) G Fa ; the integration extends over 
the part r > i? of the Cauchy hypersurface S^. At this place one inserts the bounds on \E'o, 
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and \I/2- That requires some care; the estimates hold true on the initial hypersurface Eq, 
while here one needs estimates on S^. This point is clarified later. It is useful to introduce 
dimensionless variables x = R/a and fh = m/a. 

i) First we shall consider the contribution that is due to 3\l/i. Let tq be defined by 
(r, t) G Tj-Q. The insertion of the bound given in ( |3.4| ) bounds 2mi dr7]j.-^ \ by 

QmCi I drr]^ ^^^l°^ ] . Notice that gi{r) is an increasing function, therefore if one replaces 
gii^o) by gi{r), then the integral that appears here cannot be smaller. In this way one utilizes 
the initial information (the energy inequality (|3.4|) ) in order to control the evolution. The 
integral in question can be performed explicitly, with the result 

6»C.(/;*,;^)"^5!^£l^. (8.4) 

Here, it holds 

m^(137 - 770m + ISSOm^ - 21Q0m^ + 960m^) 



Gi(x) 

30(-l + 2m)5x8 
m3(991 - 5110m + 10840m2 - SSSOm^ - 720m^ + 3360m5) ^ 



^2 



m 



105(-1 + 2m)5x7 
(2981 - 13010m + 18440m2 ^ 'jg20m^ - 41520m4 + 27360m5) 



420(-l + 2m)5x6 
m(4497 - 11370m - 21720m2 _^ I33040m3 _ 200240m^ + lOlGOOm^) 

2100(-l + 2m)5x5 
375 + 4650m - 35400m2 ^ 93200m3 - llOOOOm^ + 49376m^ 
3360(-l + 2m)5x4 
11 11 11 llln( 



+ 



I ' ~2rh+x ' 



336mx'^ 448m^x^ 448m^2; 896m'^ 

T„ / —2rh+x 

(280m^ - QAOrh^x + 560m V - 224mx^ + 35x^) ^^^^^ )■ (8.5) 
This rather long expression is quite well approximated by Gi = (H-41nx)/(16x'^) ifm/a << 



1. The integration of ( |8.4| ) along a null cone Ca is done as follows. The integral \jGi{x) 
is bounded from above, 

1/2 / \ 1/2 



/OO j / rco \ f \ 

dx7]-^^Gi{x) <\j^ dxx'^ Gi{x)\ dxx-'^r]-^\ 



(8.6) 



Numerical integration yields 
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2^ 



8.24 V ^ ' ^ +0{im/af). (8.7) 
a 

One can check that the neglected terms can give a contribution comparable to the leading 
term only at distances smaller than 6.6m. 

ii) The calculation concerning the contribution of the \l/2 function is similar. The leading 
(proportional to m°) term is dr^. |\E'2| is bounded in terms of g2{x). g2{x) is an 
increasing function, and a reasoning similar to what was made when discussing gi{x), leads 
to the conclusion that one can again use the initial energy inequality given by ( p.4|) . One 
finds that dr^ is bounded from above by 



4m2(-3 + 44m - 120m2 + 9Qm^) 

21(-1 + 2m)3a;7 

2m(-21 + 236m - 408771^ - 192m3 + 576m^) 
63(-l + 2m)3x6 

-21 + 88m + 480m2 - 1968m3 + 1760m^ 
105(-1 + 2m)3F 
34 31 31 



+ 



105x4 630ma;3 840m2x2 
31 , 311n(^^, 



840771% 1680m4 

8m(60m2 - 70mx + 21^^) ln(^^^) 



(8.8) 



105x7 

In the limit of m — the function G2{x) coincides with (—4 + 5x)/(20x^). Similarly as 



before, in order to get a term bounding \/Hm, one should integrate ^jG2{x)/{l — 2m/ x) 
along a null cone Ca- That gives 0.15, up to terms 0{m), after manipulations similar to those 
done earlier. The 0{m) correction becomes dominant when 2m/ a > 0.3. After a reasoning 
similar to that applied above in the case of one finds that the total contribution due to 
the bound on the \l/2 function is equal to 



1/^ (-4 + 5x) mJEia, 0) m\ E(a,0) 

dx ^ on 3 =2-19 ^ ' ' \ (8.9) 

V Ji 2{)x'^ a a 

17 



In summary, one obtains 



Hm < 10.43 V ^ ' ^ ^ Q^^2y 



(8.10) 



In the above analysis, in ( p.3|) , we neglected the terms proportional to 63m. They give 
corrections of the order O(m^) to the right hand side of ( ^.10 ). We checked that their 



contribution is small in the region a > 6.6m. Our final result ( I.IC ) tells us that a > 
V218m 15m is valid for a nontrivial estimate. Therefore, all higher order terms in ( ^.lU ) 
can be safely neglected. 
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